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Abstract 

We consider invariant matrix models with log-normal (asymptotic) weight. It is known 
that their eigenvalue distribution is intermediate between Wigner-Dyson and Poissonian, 
which candidates these models for describing a system intermediate between the extended 
and localized phase. We show that they have a much richer energy landscape than expected, 
with their partition functions decomposable in a large number of equilibrium configura¬ 
tions, growing exponentially with the matrix rank. Within each of these saddle points, 
eigenvalues are uncorrelated and confined by a different potential felt by each eigenvalue. 
The equilibrium positions induced by the potentials differ in different saddles. Instan- 
tons connecting the different equilibrium configurations are responsible for the correlations 
between the eigenvalues. We argue that these instantons can be linked to the SU (2) compo¬ 
nents in which the rotational symmetry can be decomposed, paving the way to understand 
the conjectured critical breaking of U{N) symmetry in these invariant models. 
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I. INTRODUCTION 


Random matrices were introduced to describe the spectrum of heavy nuclei [l| , 
but over the years it has been realized that they have the power to describe a large 
universality class of many-body strongly interacting systems, with applications rang¬ 
ing from nuclear physics to string theory, from 2-D quantum gravity to condensed 
matter physics, from statistical physics to econo-physics, from neuroscience to chaos 
theory, from number theory to integrable systems, and so Although in dif¬ 

ferent applications the matrices can have different interpretations, we can use the 
partition function and the free energy as unifying concepts 

Z = Cn j = InZ , (1) 

since, as usual, they can be used as generators for all physical observables (Cat is 
an normalization constant introduced for later convenience). Here g is a coupling 
constant, M is an x Hermitian matrix and 1U(M) is the (matrix-valued) ac¬ 
tion. One can also consider different symmetry classes for the matrices, for instance 
orthogonal or symplectic, but we will concentrate just on the Hermitian case. 

Typically, matrix models are considered in a basis invariant way, that is, the 
action in ([T]) is assumed invariant under a U{N) rotation of the matrix; 1U(M) = 
trl/(M). This choice is often motivated by physical reasons, but in any case in¬ 
variant matrix models enjoy several powerful techniques for their study, which pro¬ 
vides even non-perturbative results: from the use of orthogonal polynomials to the 
Riemann-Hilbert formulation, from the Coulomb gas approach to the mapping into 
Toda lattices, and so forth jJQ. 

The versatility of random matrices is a consequence of a large universality, cap¬ 
tured by a quadratic potential and valid for any polynomial V{x). This universality 
is signaled by the correlated eigenvalue statistics known as Wigner-Dyson (WD) 
statistics. Its main feature is level repulsion: a characteristic suppression of the 
probability of finding two eigenvalues arbitrarily close. To capture this effect, one 
considers the level spacing distribution, that is, the distribution of the distance be¬ 
tween eigenvalues. All standard, invariant matrix models present the same curve, 
hallmark of the universal WD level spacing, which is shared by a variety of strongly 
interacting systems. 
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The WD statistics for the eigenvalues has emerged as a defining signature of 
extended states for the eigenfunctions. Physically, one understands this connection 
because the eigenvalues interact only through their eigenvectors and thus, to repel 
one-another, their eigenstates need to have a finite overlap. Mathematically, this 
connection is well captured by the invariant matrix models, since the eigenstate and 
eigenvalue contributions decouple in the partition function; 
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( 2 ) 


where the Hermitian matrix has been decomposed in its eigenvalues and eigen¬ 
vectors as M = IJlAU, with A a diagonal matrix of entries \j and where 
A(A) = ~ Monde determinant which emerges as the 

Jacobian in changing integration variables from the matrix elements to the eigen¬ 
value/eigenvector representation. We remark that this Jacobian introduces an ef¬ 
fective interaction between the eigenvalues, which leads to a Coulomb Gas picture: 
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in which the eigenvalues are interpreted as one-dimensional particles, confined by 
the potential V{x) and interacting through a logarithmic repulsion. The balance 
of these two contrasting forces results in an equilibrium configuration for which the 
eigenvalues spread around the minima of V{x). 

The factorization in (|2]) has lead to the conclusion that the eigenvector distri¬ 
bution is independent from V (M) and thus uniform over the A-dimensional sphere 
spanned by the Hilbert space. This behavior characterizes invariant matrix models 
as describing extended systems, i.e. systems where in any energy window at least one 
of the eigenvectors has all finite entries. This means that, when the matrix model 
is applied to mesoscopic systems, it only describes conducting states, i.e. situations 
where the disorder is not strong enough to localize the wavefunctions. 

Recently, it was realized and shown that this expectation is violated when the 
eigenvalue distribution deviates from the WD j^. In particular, if two eigenvalues 
are separated by a gap, their eigenvectors cannot mix. Without reference to a pre¬ 
ferred basis, localization means that an eigenvector can move only within a portion 
of the whole Hilbert space and that an arbitrary, not-too-strong, perturbation can¬ 
not change this. If the eigenvalues of an invariant matrix model distribute over n 
disconnected intervals, each of which containing rrij {j = 1.. .n) eigenvalues, each 
eigenvector is localized over the mj-dimensional portion of the Hilbert space spanned 
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by the eigenvectors of the other eigenvalues in that interval. This phenomenon cor¬ 
responds to a spontaneous breaking of U{N) rotational symmetry into 11^=1 
As in any spontaneous symmetry breaking (SSB), the directions of these patches of 
Hilbert space are random and uniform (they can be induced by an external small 
symmetry breaking term), but once the system has chosen them, it is hard to change 
them and impossible in the thermodynamic limit. 

In the SSB of rotational symmetry was proven for matrix model with poten¬ 
tials able to separate the eigenvalues in disconnected patches, the so-called multi-cut 
solutions. We are going to show that there exists another class of invariant matrix 
models, where the external potential is so weak that it combines with the interpar¬ 
ticle interaction to generate an effective potential induced self-consistently by the 
eigenvalue distribution. The net effect of this process is that the partition function 
can be decomposed in the contributions of several saddle points, corresponding to 
different equilibrium configurations and thus of different effective potentials felt by 
the eigenvalues. 


II. THE WEAKLY CONFINED MATRIX MODEL 


We will consider a U{N) invariant ensemble of random matrix characterized by 
a non-polynomial confining potential j^, 11|, that is, with a potential of the form 
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where n controls the strength of the interaction. This soft confinement justifies the 
name Weakly Confined Matrix Models (WCMM) for these systems. 

The logarithmic behavior sets these potentials apart from the usual matrix mod¬ 
els, where the potentials are polynomial in the matrix eigenvalues. Most importantly, 
this difference implies that the WCMM do not belong to the WD universality class. 
In fact, the level spacing is intermediate between the WD (approached for k —)■ 0) 
and the Poissonian distribution (corresponding at the k ^ oo limit), typical of 
a completely localized system Q. Moreover, the density-density correlation func¬ 
tion of the WCMM (in the unfolded coordinates) matches the expectations for the 
Anderson model at the metal/insulator transition Q, at least in the perturbative 
regimes. 

These observations have lead over the years to the conjecture that the WCMM can 
spontaneously break rotational invariance. In it was shown that in the k ^ oo 
limit the eigenvalues crystallize on a regular lattice, so that they do not interact 
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anymore. This observation was used in j9j to argue that this eigenvalue behavior re¬ 
flects the complete localization of eigenvectors, corresponding to a breaking of U{N) 
symmetry into a consistently with what expected from the level spacing be¬ 

havior. Starting from these results, we will show what happens to the eigenvalue 
distribution at finite k and argue how this behavior reflects a critical breaking of 
U{N), consistent with the conjectured relation between WCMM and the Anderson 
Metal/Instulator Transition (MIT). 

The Riemann-Hilbert problem associated to potentials of the type (jl]) cannot be 
solved. Moreover, these potentials have an indeterminate moment problem, that is, 
the same set of orthogonal polynomials have the same set of moments with respect to 
several different measures. Nonetheless, the WCMM can be exactly solved using the 
technique of orthogonal polynomials and in this way all the spectral characteristics 
have been calculated. The WCMM were actually introduced as the first example 
of a matrix model solvable using g-deformed (Hermite or Laguerre) polynomials, 
where the deformation parameter is determined by k Q. We stress once more that, 
because of the indetermination of the moment problem, different potentials sharing 
the same asymptotic behavior (jlj) are solved by the same set of polynomials. 

For the moment, let us assume that the potential is exactly log-square and that 
the eigenvalues are confined to the positive semi-axis. We perform the change of 
variables 

A, = e-^' , (5) 


in the partition function, getting 

r ^ 

z = CnVoI [f/(A^)] / Y[dxj JJ (e'^^" - , (6) 

j=l n<m 

where Vol [U{N)] = nn=i is the volume of the unitary group. 

Working in a simplex where the eigenvalues are sorted in increasing order [xj > xi 
for j > 1), we can write 


r ^ 

z = CnVoI [U{N)] / JJdXj Y[[l- ^ 

j=l n<m 

In Q it was noted that in the k, ^ oo limit, all the exponentials in the Van der 
Monde determinant can be neglected and the eigenvalues crystallize: 


Urn Z = CNyol[U{N)]K^N\exp 
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In the large k limit, the eigenvalue distribution factorizes and the interaction be¬ 
tween the eigenvalues can be exactly accounted by an new effective potential, which 
has the net effect of localizing each eigenvalue on a perfect lattice with a quadratic 
potential of width In Q| this phenomenology was discussed in details and it was 
shown that the eigenvalues become effectively uncorrelated in the k ^ oo limit. In 
Q it was further shown that this crystallization of the eigenvalues tends to localize 
the eigenvectors as well. In fact the largest eigenvalue becomes exponentially larger 
than its closest neighbor and this means that rotating the diagonal matrix with 
a generic, uniformly distributed, unitary matrix produces a matrix with negligible 
off-diagonal elements. One can thus trace out the row and column corresponding 
to this eigenvalue and the same procedure can be applied again to the new highest 
eigenvalue. This structure was further conhrmed by a numerical analysis in which 
a set of random matrices was generated with the crystallized eigenvalue distribu¬ 
tion. These matrices were perturbed by multiplying the off-diagonal elements by 
a factor slightly less than unity. The resulting eigenvector distribution showed a 
signihcant deviation from the Porter-Thomas, indicating a localization of the eigen¬ 
vectors. These arguments were the hrst direct proof that the crystallization of the 
eigenvalues affects the distribution of eigenvectors and tends to localize them. We 
can characterize this phenomenon as a breaking of the original U{N) symmetry into 
a U^{1). 


III. FINITE K ANALYSIS 


We saw that in the k ^ oo limit, the particles realize a sort of Mott-insulator 
that pins them to a lattice. In ([7]) let us set 


Cm = 


Af(4V2-l) 


(27r^)^/^Vol [U{N)] 


and dehne Xn = 2n — 1 + tn- 


(9) 


Z = f U [1 , (10) 

j=l n<m 

If the eigenvalues are ordered, for k ^ oo, the term vanishes. It can 

be checked (and it will become more apparent as we proceed) that by removing the 
square from the Van der Monde determinant, we also remove the terms responsible 
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for the other ordering of particles. Thus, in this way can work within one simplex 
without worrying about particle ordering anymore: 


Z = 



K 

2 


( 11 ) 


For large k we can now use as an expansion parameter. In fact, we see that 
the effect of each term in the Van der Monde is to add a linear displacement to the 
quadratic potential, and thus it changes the equilibrium position of the particles. 
For n = 1, a pair of nearest neighboring eigenvalues move from 2n — 1 and 2n+ 1 to 
2n. For n = 2 next-to-nearest neighbors move from 2n — 1 and 2n + 3 respectively 
to 2n and 2n + 2, and so on. After integration, each Gaussian integral contributes 
with a factor Thus, we introduce the expansion parameter q = e~'^. 

At zeroth order, from ([8]), of course we have 




r ^ 
i=i 


( 12 ) 


At first order in q we have 


N ^ N 


m=2 j=l 
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(13) 


At the next order things start to complicate a bit, because we have three types of 
contributions; from two pairs of distinct nearest neighbors, from triplets of nearest 
neighbors and from pairs of next-to-nearest: 


Zo = 
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We see that the last two contributions amount to the same process, but with equal 
and opposite weight. 
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Additional terms can be worked out similarly, but it is easy to see that the final 
result is 


N-l 


z = n (1 - 9 ”) 


N-n 


n=l 


= 1-{N-I)q + 




(16) 


which coincides with the exact result obtained in [1^, using the technique of orthog¬ 
onal polynomials (in this case, the Stieltjes-Wigert polynomials), once the overall 
normalization is taken into account. We see that the first line of 0151) derives 
directly from the Van der Monde in (ITT]) and allows for a straightforward interpre¬ 
tation of the partition function. 


IV. THE ENERGY LANDSCAPE 

Starting from the k ^ oo configuration, in which each eigenvalue Xj has equilib¬ 
rium position at — 1 (j = 1... N), each term in the Van der Monde ads a linear 
contribution to the quadratic potential and thus shifts the positions pairwise closer 
by one unity. All the configurations reached this way are stable by construction 
(small fluctuations are confined by the potential) and thus constitutes saddle points 
of the partition function. The terms of the Van der Monde connecting the different 
configurations have a natural interpretation as an instanton. 

Every equilibrium configuration with center of mass at m is realized by the action 
of the instantons; all these configurations have the same leading weight (given by 
C]sf and growing like ei'^'^^) and differ only by subleading contributions, with a 
power of q equal to the number of steps each eigenvalue traveled to reach this 
configuration, starting from the k ^ oo one. In some of these configurations, the 
equilibrium point for more than one eigenvalue can coincide (in fact, the action of 
every allowed instanton brings all eigenvalues to oscillate around the same point m). 
Thus, there is an exponential number (in N) of possible equilibrium configurations. 
In each of them, eigenvalues are only subjected to the effective potential and do not 
interact. Hence, the eigenvalue correlation is a purely instantonic effect. In fact, 
these instantons have a fermionic nature, in that each of them contributes with a 
minus sign and only one instanton connecting two given eigenvalues is allowed. 

The partition function (lT5l) has a natural interpretation in terms of the instantons, 
since only the largest eigenvalue can be connected to the smallest by a instanton with 



weight while there are two ways to connect a pair of eigenvalues N — 2 apart 

and so on, concluding with the — 1 instantons connecting nearest neighboring 
eigenvalues. 

This newly discovered energy landscape is characteristic of a complex system and 
calls for the usual tools employed in the analysis of these systems. Intriguingly, the 
mathematical framework we just outlined shares many similarities with a recently 
discussed glassy system, of hard spheres in inhnite dimensionality 12| . 

In Q, it was argued that instantons like the ones we discussed are responsible 
for restoring the broken symmetries. Since the k —)■ cxo conhgurations correspond 
to a breaking of U{N) symmetry into U{1)^, each instanton can be interpreted as 
restoring the SU{2) component connecting the two eigenvalues. 


V. CONCLUSIONS 


We showed that the partition function of the weakly conhned matrix models can 
be decomposed into a exponential number of independent equilibrium conhgurations, 
corresponding to different equilibrium points for each eigenvalue of the matrix. The 
different conhgurations are connected by instantons, which can be interpreted as 
generators of the the symmetries in which the U{N) invariance is broken, according 
to fl. 

The WCMM is conjecture to represent a toy model for the Anderson Metal 
Insulator Transition. The structure unveiled supports this conjecture and draws a 
novel connection between localization problems and systems with complex energy 
landscapes, such as the Sherrington-Kirkpatrick model 13| . 

Moreover, the analytically continued WCMM is connected to Chern-Simons the¬ 
ories 10| and SUSY string theories, such as ABJM [l5|. In particular, the ABJM 


partition function in the localization limit corresponding to the matrix models can 
be interpreted in terms of a free Fermi gas at hnite temperatures |3|: it seems that 
this is a gas of instantons and thus the rich energy landscape should be analyses in 
this held theory formulation as well. 

We considered models with an exact log-normal potentials: more general poten¬ 
tials can be considered, but we believe that the qualitative structure of the energy 
landscape will not change. Such generalizations include potentials which can be 
written as polynomials in logarithm of the matrix. A more interesting extension 
will be the treatment to the case of eigenvalues over the whole real axis, where 
non-local correlations between eigenvalues at opposite sides of the origin have been 
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and interpreted in terms of a Lnttinger liqnid in a Hawking-Hartle 
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